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ABSTRACT 

Recent spectro-polarimetric observations of a sunspot showed the formation of 
bipolar magnetic patches in the mid penumbra and their propagation toward the 
outer penumbral boundary. The observations were interpreted as being caused 



by se a-serpent magnetic fields near the solar surface (jSainz Dalda fc Bellot Rubio 



20081 ). In this Letter, we develop a 3D radiative MHD numerical model to explain 
the sea-serpent structure and the wave-like behavior of the penumbral magnetic 
field lines. The simulations reproduce the observed behavior, suggesting that the 
sea-serpent phenomenon is a consequence of magnetoconvection in a strongly in- 
clined magnetic field. It involves several physical processes: filamentary structur- 
ization, high-speed overturning convective motions in strong, almost horizontal 
magnetic fields with partially frozen field lines, and traveling convective waves. 
The results demonstrate a correlation of the bipolar magnetic patches with high- 
speed Evershed downflows in the penumbra. This is the first time that a 3D 
numerical model of the penumbra results in downward directed magnetic fields, 
an essential ingredient of sunspot penumbrae that has eluded explanation until 
now. 
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Introduction 



It is well-known that the sunspot penumbra (the outer part of sunspots) has a very 
complicate filamentary structure and a strong non-stationary outflow. This outflow is re- 
sponsible for t he so-called Ey ershed effect, a Doppler shift of the spectral lines emerging 
from sunspots (jEvershedlll909l ) . The magnetic structure of the penumbra can be represented 
as a mixture of two magnetic fi e ld components with different inclinations and strengths 



(e.g., 


Deeenhard 


t & Wiehr 


1991 




Schmic 


t et a. 


. 1992 




Title et a 


199; 


3; 


Lites et al. 


1993; 


Stanc 


lfiek 


et al. 


1997; 


Bellot Rubio et al. 


2004: 


Sanchez Almeida 


2005; 


Borrero et al. 


2005; 


Beck 


2008 


). However, the radial Evershed flow (in particular, the high-speed 'Evershed 



coluds') is associated w i th the more strongly incli ned, almost horizontal field (e.g. iTitle et al. 
19931 : IShine et al.lll994j ; bellot Rubio et alJl2003h . 



Recently, significa nt progress in our understanding of the Evershed effect was made by 
numerical simulations (IHeinemann et al.l 120071 : IScharmer et al.l 12008c iRempel et al.ll2009al lb; 
Kitiashvili et al.ll2009al ). These studies suggested that the Evershed flow is a consequence 
of overturning magnetoconvection in the presence of inclined magnetic fields, and that the 
driving mechanism is associated with tr aveling convective waves that propagate in t he di- 



rection of the magnetic field inclination ( jHurlburt et al. 1 120001 iKitiashvili et al. 



issue is not resolved and other interpretations are possible (see ISchlichenmaier 



2009ah . The 
2009h . Thus, 



it is important to confront the simulations with as many observations as possible. 

In th is Letter, we examine the idea th at the sea-serpent penumbral field lines de- 
tected by ISainz Dalda &; Bellot Rubio! (120081 ) in high-resolution Hinode measurements are 
related to the same mechanism of overturning convection and traveling convective waves in 
a strong inclined magne tic field. Our analysis is based on the radiative MHD simulations of 
Kitiashvili et al.l J2009ah . 



Moving bipolar magnetic patches in the penumbra and the sea-serpent 

model 



Analyses of visible and near-infrared spectro-polarimetric observations from the ground 
revealed that in the mid and outer penumbra the Evershe d flows are directed downward 
and have magnetic polarities opposite to that of the spot (IWestendorp Plaza et al.l 119971 . 
see Bellot Rubio 2009 for a detailed descri ption). The spectro-polarimeter aboard Hin- 
ode (IKosugi et al.l 120071 : iTsuneta et al.l 120081 ) made it possible to identify isolated magnetic 
patches of opposite polarity associated with strong, even supersonic, Ever shed downflows 
( llchimoto et al.ll2007l ; iBellot Rubioll2009l ; ISanchez Almeida fc Ichimotoll2009l ) . This provided 
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a beautiful confirmation of the results o btaine d from ground-based observations at lower 
resolution. ISainz Dalda fc Bellot Rubiol ( 120081 ) then discovered that the opposite polarity 
patches move radially outward, often accompanied by anoth er patch having the polarity of 



the sp ot a nd being located further away from the umbra. ISainz Dalda fc Martinez Pillet 



(120051 ) and iRavindra I tood ) also detected t he evolution of rever sed-polarity patches in the 
penumbra using SOHO/MDI observations (IScherrer et al.lll995l ). but their bipolar nature 
could not be established. Taken together, these results provide an indication that the mov- 
ing magnetic features observed in the penumbra and the Evershed flow belong to the same 
physical mechanism. 

The magnetic patches move between the more vertical filaments of the penumbra with 
a typical velocity of 0.3 — 1 km/s, their l ength is 2-3 arcsec, the mean widt h is ~ 1.5 arcsec, 
and they have a lifetime of 0.5 - 7 hours (ISainz Dalda fc Bellot Rubioll2008l ). Sometimes, the 
bipolar pairs appear one behind another along the same filament. This can be the signature 
of a wave-like behavior, similar t o the wave behavior of magnetoconvection in the penumbra 
found in numerical simulations ( iKitiashvili et al.ll2009a! Jbl). 



Based on the observational data, ISainz Dalda fc Bellot Rubiol (120081 ) proposed a sea- 
serpent model which explains the existence of moving bipolar structures as wavering magnetic 
field lines. In this model, almost horizontal field lines return to the solar interior and then 
come back to the photosphere, thus forming bipolar patches. This scenario is similar to 



that discussed by Schlichenmaier 
similar to the model prop osed by 



2002) in the context of the moving tube model. It is also 



the sunspot moat (MMFs; ISheeley 



Harvey fc Harveyl ( 119731 ) for moving magnetic features in 



19691 ). Indeed, it has been suggested that MMFs represent 



Ravindra 2006; Sainz Dalda & Bellot Rubio 2008; Kubo et al. 



ra ( 


Sainz Dalda & Martinez Pillet 


2005 


Kubo et al. 


20081). 



3. Simulations of magnetic structure and dynamics of sunspot penumbra 



For numerical sim ulations of sunspot penumbra conditions we use the 3D radiative 
MHD code "SolarBox" (IJacoutot et al.ll2008al ; IJacoutot et al. Il2008bl ). This code takes into 
account all essential physics and includes sub-grid scale turbulence modeling based on the 
large-eddy simulation (LES) approach. A dynamic Smagorinsky turbulence model provides 
the best agreement with observations in terms of the acoustic oscillation power. However, 
in these simulations we use a computationally more efficient hyperviscosity model, because 
it shows results qualitatively very similar to the dynamic model. The code uses a real-gas 
equation of state, takes into account ionization and excitation of all abunda nt species in the 
LTE approximation, and includes radiative transfer and magnetic effects (jRipoli fc Wray 
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20031 : IJacoutot et al.ll2008al : iKitiashvili et al.l 12009 



We simulate the behavior of solar magnetoconvection for an initial background magnetic 
field of various strengths, from 600 to 2000 G, but with the same inclination of 85° to 
the vertical axis. Our previous re sults have shown that these condition provide a very 
good model for the Evershed effect (jKitiashvili et al.ll2009ai ). The computational domain is 
6.4 x 6.4 x 6 Mm 3 and the grid resolution is 50 km. The results were verified using a higher 
resolution of 25 km and a larger horizontal size of 25 Mm. The lateral boundary conditions 
are periodic, and we keep the mean initial inclination of the field in the whole box domain 
by setting up top and bottom boundary conditions that conserve the total magnetic flux, 
but locally the field strength and inclination can freely change. 

The magnetic field strongly influences convective motions. For instance, in the case of 



(Stein & Nordlundl 


2001) 


(Jacoutot et al. 


2008b|). 



A strong inclined magnetic field deforms the granules, making them 



elongated along the field lines. When the field becomes almost horizontal the convective 
cells form a filamentary structure, and hi gh-speed flows (reaching 5-6 km/s) appear in the 
direction of the magnetic field inclination (jKitiashvili et al.ll2009af ). These flows share many 
similarities with the observed Evershed flows, including their nearly-supersonic speeds. 

For illustration of our results we use a simulation run with initial magnetic field of 
B = 1200 G and inclination angle of 85° pointing in the direction of the horizontal x-axis. 
Figure [1] shows a time sequence of 2D snapshots of the horizontal and vertical components 
of the magnetic field, B X ,B Z , and velocity, v x ,v z , with 8-min cadence for five time moments: 
0, 8, 16, 24 and 32 mins. The overplotted curves are projections of the magnetic field lines. 
The black and white areas seen in the vertical magnetic field (Fig. Hk) correspond to patches 
of negative and positive polarities. The horizontal magnetic field component (Fig. [Tb) shows 
variations of the magnetic filamentary structure, which is a result of magnetoconvection in 
the inclined field. The dark areas correspond to "magnetic g aps" , areas with a weaker (but 



not z ero) magnetic field, that correspond to observations (jSevernyil Il965t llchimoto et al. 



20071 ). It is clear that the areas of weaker horizontal magnetic field correspond to the tran- 
sitions between negative and positive polarities of the vertical field. In these transition 
regions, strong downward flows are observed (Fig. [Tb). The transition regions are also char- 
acterized by strong ~ 4 — 6 km/s horizontal velocities (Fig. [Tbl). The strongest horizontal 
flow in the direction of the magnetic field inclination is associated with the negative polari- 
ties. The correlation between downflows and negative polarity patches is in agreement with 
the observations, which show that the d ownward penumbral motions occur along magnetic 
fields returning back to the solar surface ( jWestendorp Plaza et al.lll997l ; llchimoto et al.ll2007l ; 
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Bellot RubiolBOOa ISanchez Almeida & Ichimotol 120091 ). 



Figure [2] displays the 3D structure of magnetic field lines (red curves), temperature 
variations, and the velocity field (arrows). The top horizontal slice corresponds to the solar 
surface. The results demonstrate a strong coupling between the fluid flow and magnetic field. 
In particular, the upflows correspond to magnetic field lines rising above the surface (positive 
polarity) while the downward flows correlate with the negative polarity. This corresponds 
very well to the observed properties of the penumbra and the sea-serpent model. 



4. Discussion and Conclusion 

In this Letter, we have used the results of numerical simulations of magnetoconvection in 
strong inclined magnetic field to interpret polarimetric observations of a sunspot penumbra. 
The results reproduce the moving bipolar magnetic elements observed in high-resolution 
SOHO/MDI and H inode/SOT data and also their pro perties, supporting the sea-serpent 



model proposed by ISainz Dalda fc Bellot Rubiol (120081 ). The simulations explain the sea- 



serpent structure and dynamics of the penumbral field as a consequence of solar magneto- 
convection in a highly inclined, strong magnetic field, which forms filamentary structures 
and has properties of traveling convective wave. 

The physical picture schematically illustrated in Figure [3] is the following. Convective 
cells in sunspot penumbrae are deformed under the action of the incline d magnetic field. 



formin g filamentary structures and producing high-speed Evershed flows (jKitiashvili et al. 



2009al ). The magnetic field lines are stretched by the downward flows and dragged under 
the surface. The points where the magnetic field lines cross the solar surface are observed 
as magnetic patches of positive and negative polarities. Note that the negative patch is 
closer to the umbra, in agreement with the observations. The convective cells move in the 
direction of the magnetic field inclination because of the traveling convective wave behavior. 
Therefore, the bipolar magnetic patches also move in the same direction. 

Thus, the numerical simulations connect the sea-serpent structure of the moving bipolar 
magnetic pathes observed in the penumbra with the process of overturning magnetoconvec- 
tion, traveling convective waves, and the Evershed flow. 



This work was partly supported by NASA, the Center for Turbulence Research (Stanford 
University), NORDITA, AlbaNova Univ. Center (Stockholm), and the Spanish Ministerio 
de Ciencia e Innovation through project AYA2009-14105-C06-06 and by Junta de Andalucia 
through project P07-TEP-2687. 



-6 - 



REFERENCES 

Beck, C. 2008, A&A, 480, 825. 

Bellot Rubio, L.R. 2009, in Magnetic Coupling between the Interior and the Atmosphere of 
the Sun, ed. S.S. Hasan & R.J. Rutten, Astrophysics and Space Science Proceedings, 
(Berlin: Springer), in press, (larXiv :0903.3619). 

Bellot Rubio, L. R., Balthasar H., Collados, M. & Schlichenmaier, R. 2003, A&A, 403, L47. 

Bellot Rubio, L. R., Balthasar, H., & Collados, M. 2004, A&A, 427, 319. 

Borrero, J. M., Lagg, A., Solanki, S. K., & Collados, M. 2005, A&A, 436, 333. 

Degenhardt, D. & Wiehr, E. 1991, A&A, 252, 821. 

Evershed, J. 1909, MNRAS, 69, 454. 

Harvey, K. & Harvey, J. 1973, Sol. Phys., 28, 61. 

Heinemann, T., Nordlund, A., Scharmer, G. B., & Spruit, H. C. 2007, ApJ, 669, 1390. 

Hurlburt, N. E., Matthews, P. C, & Rucklidge, A. M. 2000, Sol. Phys., 192, 109. 

Ichimoto, K., Shine, R. A., Lites, B., Kubo, M., Shimizu, T., et al. 2007, PASJ, 59, S593. 

Jacoutot, L., Kosovichev, A. C, Wray, A. A. & Mansour, N. N. 2008a, ApJ, 682, 1386. 

Jacoutot, L., Kosovichev, A. C, Wray, A. A. & Mansour, N. N. 2008b, ApJL, 684, L51. 

Kitiashvili, I. N., Kosovichev, A. C, Wray, A. A. & Mansour, N. N. 2009a, ApJL, 700, 178. 

Kitiashvili, I. N., Jacoutot, L., Kosovichev, A. C, Wray, A. A. & Mansour, N. N. 2009b, 
Proc. of Int. Conf. Stelar Pulsation: challenges for theory and observation. AIP Conf. 
Proa, 1170, 569. 

Kitiashvili, I. N., Kosovichev, A. C, Wray, A. A. & Mansour, N. N. 2009c, Proc. of ParCFD 
Conf., NASA Ames, 424. 

Kosugi, T., Matsuzaki, K., Sakao, T., Shimizu, T., Sone, Y., et al. 2007, Sol. Phys., 243, 3. 

Kubo, M., et al. 2008, ApJ, 681, 1677 

Lites, B. W., Elmore, D. F., Seagraves, P., & Skumanich, A. P. 1993, ApJ, 418, 928 
Rempel, M., Schussler, M., & Kndlker, M. 2009, ApJ, 691, 640. 



-7- 



Rempel, M., Schiissler, M., Cameron, R. H., & Knolker, M. 2009, Science, 325, 171. 

Ripoli, J.-F. & Wray, A. A. 2003, in Annual Research Briefs-2003, Center for Turbulence 
Research, Stanford University, 3. 

Ravindra B. 2006, Sol. Phys., 237, 297. 

Sainz Dalda, A. & Bellot Rubio, L. R. 2008, A&A, 481, L21. 

Sainz Dalda, A. & Martinez Pillet, V. 2005, ApJ, 632, 1176. 

Sanchez Almeida, J. 2005, ApJ, 622, 1292. 

Sanchez Almeida, J., & Ichimoto, K. 2009, A&A, 508, 963. 

Scharmer, G. B., Nordlund, A., & Heinemann, T. 2008, ApJ, 677, L149. 

Scherrer, P. H., et al. 1995, Sol. Phys., 162, 129. 

Schlichenmaier, R. 2009, Space Science Reviews, 144, 213 

Schmidt, W., Hofmann, A., Balthasar, H., Tarbell, T. D. & Frank, Z. A. 1992, A&A, 264, 
L27. 

Stanchfield, D. C. H., II, Thomas, J. H., & Lites, B. W. 1997, ApJ, 477, 485. 
Sheeley, Jr. 1969, Solar Physics, 9, 347. 

Shine, R. A., Title, A. M., Tarbell, T. D., Smith, K., Frank, Z. A., & Scharmer, G. 1994, 
ApJ, 430, 413. 

Stein, R. F. & Nordlund, A. 2001, ApJ, 546, 585. 

Severnyi, A. B. 1965, Soviet Astronomy, 9, 171. 

Title, A. M., Frank, Z. A., Shine, R. A., Tarbell, T. D., Topka, K. P., Scharmer, G. & 
Schmidt, W. 1993, ApJ, 403, 780. 

Tsuneta, S., Ichimoto, K., Katsukawa, Y., Nagata, S., Otsubo, M., et al. 2008, Sol. Phys., 
249, 167. 

Westendorp Plaza, C, del Toro Iniesta, J. C, Ruiz Cobo, B., Martinez Pillet, V., Lites, 
B. W., & Skumanich, A. 1997, Nature, 389, 47. 



This preprint was prepared with the A AS IATgX macros v5.2. 



-8 - 




Fig. 1. — Snapshots of the vertical (a) and horizontal (b) components of magnetic field, 
and the vertical (c) and horizontal (d) components of velocity at the photospheric level for 
different moments with a cadence of 8 min for the simulations for the initial Bq = 1200G 
background magnetic field with the 85° inclination angle from the vertical in the z-direction. 
The overlayed curves are projections of magnetic field lines. 
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Fig. 2. — 3D rendering of temperature variations, magnetic field lines (red curves) and flow 
field (arrows). 



magnetic field line 




convective cells 



Fig. 3. — Schematic illustration of the sea-serpent model of sunspot penumbra and the 
moving bipolar magnetic patches. The inclined elliptical curves illustrate streamlines of 
magnetoconvection in highly inclined magnetic field. The thick solid curve illustrates the 
sea-serpent behavior of magnetic field lines. Dashed line indicates the solar surface, where 
the magnetic field is observed. 



